Non-Markovian dynamics of a single electron spin coupled to a nuclear spin bath 
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We apply the time-convolutionless (TCL) projection operator technique to the model of a central 
spin which is coupled to a spin bath via nonuniform Heisenberg interaction. The second-order 
results of the TCL method for the coherences and populations of the central spin are determined 
analytically and compared with numerical simulations of the full von Neumann equation of the total 
system. The TCL approach is found to yield an excellent approximation in the strong field regime 
for the description of both the short-time dynamics and the long time behavior. 
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I. INTRODUCTION 

Projection operator techniques^ are widely used in 
studies of the dynamical behavior of complex open quan- 
tum systems featuring non-Markovian relaxation and de- 
coherence phenomena^. The most prominent variant of 
these techniques is the Nakajima-Zwanzig (NZ) projec- 
tion operator method which leads to an integrodifferen- 
tial equation for the reduced density matrix of the open 
system containing a certain memory kernel^. An alter- 
native and technically much simpler scheme is the time- 
convolutionless (TCL) projection operator technique in 
which one obtains a first-order differential equation for 
the reduced density matrij*^. The advantage of the TCL 
approach consists in the fact that it yields an equation 
of motion for the relevant degrees of freedom which is 
local in time and which is therefore often much easier 
to deal with than the NZ master equation. In fact, this 
method has been applied to many physical systems show- 
ing strong non-Markovian effects (see, e. g., Refs.HQS)- 

In the present paper we apply the TCL projection op- 
erator technique to the model of a central spin interacting 
with a bath of N spins defined by the Hamiltonian {11=1) 



H 



N 



-(73 



k=l 



(1) 



The Pauli operators cr and cr'^ act on the Hilbert spaces 
of the central spin, which is regarded as the open quan- 
tum system, and of the k-th bath spin, respectively. The 
strength of the spin-bath coupling is given by the con- 
stants ak- Moreover we have included an external mag- 
netic field that acts on the central spin and leads to the 
Zeeman splitting loq- 

The model given by the Hamiltonian ^ may be used 
to describe for example a single, localized electron spin 
coupled to a bath of nuclear spins in a quantum dot 
through contact hyperfine interaction^.. It features many 
interesting phenomena such as non-exponential behavior 
of correlations and coherences and strong non-Markovian 
effects. A detailed treatment of the model within the 



NZ projection operator technique has been carried out 
m Ref. |10|, and non-perturbative solutions for polarized 
initial conditions have been constructed in Ref. [111. Re- 
cently, a detailed analytical and numerical study of the 
exact Bethe ansatz solutionis of the model has been car- 
ried ouijii. Moreover, several efficient numerical algo- 
rithms have been proposed that are based, e. g., on the 
spin-coherent-state reprcsentationi^ or on the Chebyshev 
expansion for the full propagator—. 



The TCL projection operator method has been ap- 
plied to various spin bath models for which a com- 
pact analytical solution is available, such as central spin 
models with Heisenberg XY— and with full Heisenberg 
intcractiortii for uniform couplings, and central spin 
models with nonuniform Ising intcractioni^. The pur- 
pose of the present paper is a detailed investigation of 
the performance of the TCL technique for the nontrivial 
model given by Eq. ^ with nonuniform couplings. To 
this end, we will compare the results for the populations 
and the coherences of the central spin obtained from the 
TCL approach with numerical simulations of the full von 
Neumann equation of the model. It will be demonstrated 
that the method provides an efficient scheme which is 
applicable in the perturbative regime of weak couplings, 
even for long interaction times. 



The paper is organized as follows. Section [III contains 
a brief account of the TCL projection operator technique 
and its application to the model given by the Hamilto- 
nian ([T]), as well as the derivation of the master equation 
governing the dynamics of the reduced density matrix of 
the central spin. We compare in Sec. Illll the solutions of 
this master equation with numerical simulations of the 
von Neumann equation corresponding to the Hamilto- 
nian In Sec. IIVI we discuss the performance of an 
alternative TCL approach that is based on a modified 
interaction picture and leads to a simplified master equa- 
tion. Finally, we draw our conclusions in Sec. IVl 
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II. TCL MASTER EQUATION 
A. Interaction picture 

It is convenient to write the Hamiltonian ((T|) as H = 
Hq + Hi, where 

Ho = + 2^3X3 (2) 

represents the unperturbed part, and 

Hi = 2{a+K_+a_K+) (3) 

is the interaction HamiltonianiS. In Eq. ([3]) (j± are the 
raising and lowering operators of the central spin whereas 

^ N N 

K3 = 2 ^ "'^'^3 , = ^ a^ai . (4) 

k=l k=l 

In the interaction picture defined by Hq the interaction 
Hamiltonian becomes 

Hiit)^a+B_it) + a.B+it), (5) 

where 

B±{t) = 2eT''^''*e=F2'■^^3t^^gT2^A:3t, (g) 

Thus the dynamics of the total system's density matrix 
p{t) is governed by the von Neumann equation 

j^pit) = -,[Hi{t),pit)]^Cit)p{t), (7) 

where C{t) denotes the Liouville superoperator corre- 
sponding to the interaction Hamiltonian Hi(t). 

B. TCL projection operator approach 

The starting point of the projection operator technique 
is the introduction of a suitable projection superoperator 
V. This is a positive and trace preserving linear map 
that acts on the operators of the total system with the 
property of a projection operator, i. e. = V. The 
superoperator V is used to project any state p of the total 
system onto its relevant part Vp, expressing formally the 
elimination of the irrelevant degrees of freedom from the 
full dynamical description of the underlying mode]3. 

Projection operator techniques are used to derive a 
closed equation of motion for the relevant part Vp. A spe- 
cial variant of these techniques is the time-convolutionless 
(TCL) projection operator method. Given an initial state 
p(0) satisfying Vp{0) = p(0), this technique leads to 
a time-local first-order master equation for the relevant 
part of the form 

^Vp{t) = IC{t)Vp{t). (8) 



Here, /C(t) is a certain superoperator, representing the 
explicitly time-dependent generator of the quantum mas- 
ter equation for Vp. We note that, like the corresponding 
NZ equation, the TCL master equation ([5]) describes all 
non-Markovian effects although it is local in time. 

In practical applications the TCL generator /C(t) is 
mostly obtained from a perturbation expansion with re- 
spect to the strength of the interaction Hamiltonian, 

/C(i) -/Ci(i) + /C2(t) + ... (9) 

The first-order contribution is given by 

ICi{t)^VC{t)V, (10) 

while the second-order term takes the form 

/C2(t)=/ dti[VC{t)L{ti)V -VL{t)VC{ti)V]. (11) 
Jo 

We remark that this expansion corresponds to an expan- 
sion in terms of the ordered cumulants of the Liouville 
operator— C{t). 

In the present paper we restrict ourselves to the second 
order and employ the following projection superoperator, 

Pp = ^TrB{H™p}®-^H™. (12) 

Here, Tr^ denotes the partial trace over the spin bath and 
the H,„ are ordinary projection operators acting in the 
conventional sense on the Hilbert space of the spin bath. 
They project onto the eigenspaces of the 3-component of 
the bath angular momentum, 

1 ^ 

h = ^Y.^l (13) 
fc=i 

corresponding to the eigenvalues m — — ^ 1 ■ • ■ j y ■ "^^^ 
quantity 

7V„ ^TrsH™ = (^^^^^^ (14) 

represents the degree of degeneracy of the eigenvalue m 
of J3. Explicitly, we have 

H„i = ^m^ , . . . ,m^){m^ ,w? , . . . \, 

m^—m 

(15) 

where m'' = ±^ denotes the eigenvalue of the fc-th bath 
spin operator ^(t|. Obviously, the projection operators 
YLm fulfill the relations, 

IlmHm/ = (5mm'H,„/, ^ H^ = /. (16) 

m 

With the help of Eq. it is easy to verify that the 
projection superoperator (|12p is indeed a completely pos- 
itive and trace-preserving map that satisfiesi^ = V. 
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It projects a given state p onto a separable quantum 
state Vp which describes classical correlations between 
the (unnormalized) system states 



Pm{t) = TYB{^mp{t)} 



(17) 



and the bath states n,„/iVm- The latter represent states 
of maximal entropy under the constraint of a given value 
m for the total angular momentum. Finally, the reduced 
density matrix of the central spin is given by 



ps{t) = Tri3 p{t) = ^ p,n{t). 



(18) 



Thus, the dynamics of the central spin is determined by 
the dynamical variables pm{t), m = — ■ ■ ■ , y- 

As mentioned already the TCL master equation ([8]) 
presupposes that the total system's initial state p{0) ful- 
fills the condition 



^P(0) = P(0). 



(19) 



If this condition is not satisfied one has to add a certain 
inhomogcneity to the right-hand side of the TCL mas- 
ter equation which involves the initial conditions through 
the complementary projection Qp{0) = {I ~ V)p{0). 
In the standard applications of the projection operator 
techniques one employs a projection superoperator that 
projects onto an imcorrelated tensor product state. Con- 
dition (|19|1 then holds, of course, if and only if the initial 
state is an uncorrelated state of the form p(0) = ps®PBi 
where ps is a state of the system and is a state of the 
bath. However, the projection given by Eq. ([T^ belongs 
to the class of correlated projection superoperators^ii^ 
which projects a given state p onto a correlated system- 
bath state. For this correlated projection the condition 
(fT9|) is satisfied if and only if the initial state is of the 
form 



p(o) = ^p„(o)®-^n. 



(20) 



Hence, the initial state may contain certain statistical 
correlations. A great advantage of the correlated pro- 
jection operator technique is therefore given by the fact 
that it allows the treatment of correlated initial states by 
means of a homogeneous TCL master equation. 



C. Deriving the master equation 

For the interaction Hamiltonian ([5]) the projection op- 
erator satisfies VC{t)V = 0. Hence, to second order 
the TCL master equation dH) reduces to 



^Vp{t) = dtiVC{t)C{h)Vp{t). 



(21) 



Taking into account the definition of V as given by 
Eq. and exploiting the properties (fTB|) . it is possi- 
ble to convince oneself that Eq. (PT|) is equivalent to the 



following system of coupled differential equations in the 
dynamical variables Pra{i)i 



X Trn <^ H 



Hi{t), 



Evaluating the double commutator we finally get the 
master equation 

^Pm(i) = dT^[g,n+i{T) + gl,+i{T)]a+p,n+i{t)a_ 

+ [/m-l(T) + f*,_i{T)]a^Pm~i{t)c7+ 
-fm{T)(J+a^Pra{t) - (t)/?™ (t)(T+(T_ 

-g,n{T)(T-(J+p^{t) - 5;„(r)p™(i)cr_cr+|. 

(22) 

The correlation functions fm{T) and g,„(T) are defined 
by 



/,„(r) = (i?_(t)i3+(ti)),„, 
with T = t — ti and 

= -r^TrsjOn™}. 

Exploiting Eq. ([5]) wc find 

/,„(t) = 4^a2(a^a^e'("°+4^^+2"'^-)^ 

„j(-a;o-4A'3+2Qfc)T 



,(r) = ^Y.^lUat 



(23) 
(24) 



(25) 



(26) 



.(27) 



III. COMPARISON WITH NUMERICAL 
SIMULATIONS 



In this section we compare the dynamics of the den- 
sity matrix psit) of the central spin obtained by solving 
the master equation ([^ with the one deduced directly 
from the von Neumann equation ([7|. In particular we 
will consider the dynamics of both the coherences and 
the populations of the central spin starting from a fixed 
initial condition. To this end, we have carried out numer- 
ical simulations of the full von Neumann equation with 
mixed initial states for systems with up to A^ = 10 bath 
spins. Of course, the number of bath spins for which a 
direct numerical simulation of the von Neumann equa- 
tion is possible is limited by the exponential increase of 
the dimension of the underlying Hilbert space. We note 
that the dimension of the total state space (the space of 
density matrices) is given hy D ~ 2^^+^ — 1, which yields 
D « 4 • 10^ for AT = 10. 
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In the following we assume that the hyperfine coupling 
constants are given by 



Uk = ao exp 



^/d' 



(28) 



where fco = N/2. We choose n/d = 2 corresponding 
to a Gaussian electronic wave function [n ~ 2) in one 
dimension (d = 1) which is of relevance in the context 
of a quantum do t^°'^'^ . We denote by Ai the mean of 
the coupling constants and by A2 the respective root 
mean square, 



"fe, 



A. = 




(29) 



The initial state of the total system is taken to be 
p{Q) = ps(0) ® Pb{0), where pb(0) = 2^^I represents 
an unpolarized infinite temperature state (/ denotes the 
unit matrix of the spin bath). With this initial state 
we have Pm(0) = N,nPs{0). We emphasize that the 
present technique also allows the treatment of polarized 
and of correlated initial states (see Sec. Ill B) ). 



Coherences 



must be multiplied by the factor (e Hence, the 

coherence C{t) in the rotating frame of the central spin 
is found to be 



(34) 
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FIG. 1: (Color online) Real and imaginary part of the co- 
herence of the central spin for ao/^^o = 0.01 and A'' = 10 
bath spins. Blue dots: Numerical simulation. Red line: TCL 
approximation according to Eq. (|34|) . 



The coherence of the central spin is defined by C (t) = 
{+\ps{t)\—), where |±) denote the eigenstates of 0-3. Ac- 
cording to Eq. we have 

c(i) = Ean(<) = E(+K(OI-)- (30) 

7n rn 

Starting from the master equation (|22p we have 

jCm{t) ^ - dT[f^{T) + g*^{T)]C^{t) (31) 

with the obvious solution 

^^(O =C,„(0)e-^""W, (32) 

where 

A™^i) = /* dh r dh [fra{t2) + .9,;(i2)] • (33) 

Let's now remember that the master equation (22) has 
been written in the interaction picture with respect to the 
Hamiltonian As usual, we will represent our results 
in the interaction picture with respect to the free Hamil- 
tonian ^(T3, that is in the rotating frame of the central 
spin. In order to do this we have to use the replacement 

E /'™(^) « ^ - E e-^"^^^^V™(0 ® ^e^-3^3t. 

m m 

It is immediate to observe that under this transformation 
the populations remain unchanged, while the coherences 



In Fig. [T] we compare the result of the TCL approxi- 
mation given by Eq. (|34[) with numerical simulations of 
the von Neumann equation ([7]) supposing that the cen- 
tral system is initially prepared in the superposition state 
7f (1+) + h))- The correlation functions (gHl) and ((77)) . 
and the mean value (e""''^^*') have been calculated nu- 
merically. As it is evident from the Figures the agree- 
ment between the TCL result and the numerical solution 
is excellent in the perturbation regime, even for long in- 
tegration times. 

B. Populations 

The populations P±{t) ~ {zt\ps{t)\zk) of the central 
spin are given by 

P±it) = E^™ W = J2(^\Prn{t)\±). (35) 
m m 

The master equation ((22)) leads to a system of coupled 
equations, 

d /■* 

-[/™(r)+/;„(r)]P+(0, (36) 
J^Pmit) = I rfr[/„,_i(r)+/;„_i(r)]P+_i(0 

~[9^n{r)+g*^{T)]P~{t). (37) 
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To solve these equations we employ the relation 
d 



dt 



[Pr^^{t) + P-,+ ,it)] =0, 



(38) 



which expresses the conservation of the 3-component of 
the total spin angular momentum. Using the initial con- 
dition P+{0) = 1 we obtain 



1+ / d<ie^™^(*^)Mrn(tl) 
"'0 



(39) 



where 



A^°P(t) = 2Re f dh f dh [9rn+l{t2) + fmit2)] , (40) 

Jo Jo 



and 



fi,nit)=2Re drgm+iir). 
Jo 



(41) 




FIG. 2: (Color online) Population of the central spin for 
ao/uJo = 0.01 and A'^ = 10 bath spins. Blue dots: Numer- 
ical simulation. Red line: TCL approximation according to 
Eq. dMll- 



The comparison of the TCL result p9|) with the nu- 
merical simulation is shown in Fig. [21 We again observe 
a very good agreement of the TCL approximation with 
the exact dynamics for short and also for long interac- 
tion times. Due to the exponential increase of the nu- 
merical effort with increasing N we can treat only a rel- 
atively small number of bath spins. Thus we have an 
open quantum system, the central spin, which is coupled 
to a relatively small environment, a finite system of bath 
spins. For this reason, the conventional techniques used 
in the theory of open systems to derive a Markovian mas- 
ter equation are not applicable because they are usually 
based on an effectively infinite environment with a con- 
tinuum of bath modes. However, the TCL approximation 
scheme developed here does not require that N be large. 
The master equation is therefore valid also for a very 



small number of bath spins. We illustrate this point in 
Figs. [3] and m which show the dynamics of the coherences 
and the populations for = 6 bath spins. 
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FIG. 3: (Color online) Real and imaginary part of the coher- 
ence of the central spin for ao/uJo ~ 0.01 and N — 6 bath 
spins. Blue dots: Numerical simulation. Red line: TCL ap- 
proximation according to Eq. p4|l . 
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FIG. 4: (Color online) Population of the central spin for 
ao/uJo — 0.01 and = 6 bath spins. Blue dots: Numeri- 
cal simulation. Red line: TCL approximation according to 
Eq. ([39)1. 



The master equation ((22)) has been obtained from 
the second order of the TCL perturbation expansion. 
Of course, for much stronger system-bath couplings the 
second-order result fails, indicating the relevance of cu- 
mulants of higher-order. To illustrate this point we have 
increased ao by a factor of 10. The result for the co- 
herences and the populations is depicted in Figs. [5] and 
El respectively. We observe that the short-time behavior 
is still correctly reproduced by the second-order, while 
there are large deviations for longer times. We conclude 
from our numerical simulations that the second order of 
the TCL scheme yields a good agreement with the exact 
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dynamics for couplings up to the order of ao/wo ~ 10~^. 
Note however that the decay of the coherence C{t) and 
the corresponding decoherence time are very well repro- 
duced even for much larger couplings as can be seen from 
Figs. [3 



0.4 \ 



and 



FIG. 5: (Color online) Real and imaginary part of the coher- 
ence of the central spin for ao/wo = 0.1 and A'' = 10 bath 
spins. Blue dots: Numerical simulation. Red line: TCL ap- 
proximation according to Eq. (|34p. 




FIG. 6: (Color online) Population of the central spin for 
ao/uJo — 0.1 and TV = 10 bath spins. Blue dots: Numeri- 
cal simulation. Red line: TCL approximation according to 
Eq. dSll). 



Hi = 2^3 (Ks - A2J3) + 2 [a+K^ + a^K^ 



(43) 



represents the interaction Hamiltonian. By contrast to 
the interaction picture of Sec. [TTl here the diagonal term 
20-3/13 of the Hamiltonian ([T]) is not completely removed 
from the interaction, but we subtract only the term 
2^2(73 J3 involving an effective coupling constant A2 that 
is equal to the root mean square of the ak [sec Eqs. (fT3)) 
and ([29]) ]. Hence, the interaction picture Hamiltonian 
now takes the form 

Hi{t) = 2fT3 {K3 - A2J3) + <y+B-it) + a^B+{t), (44) 

where 

B±{t) = 2eT*"«*eT2»^=-'3t^^gT2^A2J3t^ (45) 

The corresponding Liouville operator will again be de- 
noted by C{t). 

An important point of the new interaction picture is 
that for the Hamiltonian and the projection 
the first order term VC^t)!^ does not vanish. Hence, 
one has to use the full expressions (jTU]) and (fTTI) for the 
second-order TCL generator. However, the great advan- 
tage of the present procedure is the fact that the cor- 
relation functions ((23)) and p4|) take on a very simple 
form. 



/™(r) = i3+(m)e'"+M^ 
ff™(T) = i?_(m)e^"-(™)^ 



(46) 
(47) 



where 



and 



VL±{m) = ±uoq + AA2 ±m + - 



N 

B±{m) ^AAi[—Tm 



With the help of these expressions we find the master 
equation 



IV. MODIFIED INTERACTION PICTURE 

A certain disadvantage of the perturbation scheme 
used in Sec. [ill consists in the fact that the time in- 
tegrals over the correlations functions in Eqs. ([55)) and 
([39)) - ([4T)) are generally difficult to calculate and have to 
be determined numerically. To avoid the appearance of 
these expressions and to obtain a simpler approximation 
scheme we employ a modification of the interaction pic- 
ture Hamiltonian. To this end, we write the Hamiltonian 
^ again as H = Hq + Hj, where now the unperturbed 
part is given by 



UJQ 



0-3 -I- 2^2(73 J3, 



(42) 



dt 



p.,n{t) = 2im{A2- Al)[a3,Pr,^{t)] 



^2 _ 4j^2 



t 



1 



(Al-Al)t[a3,[a3,p^{m 



dT 



|_B_(m + 1)2 cos[f2_|_(TO)r]o'+p,„+i(t)cr_ 
+i3+(m — l)2cos[r2_(m)T]CT_p„i_i(t)(T+ 
-B+(m)e'"+(™)V+a_p,„(t) 



— i3+(m)e 



Prn{t)i7+a- 



-B_(m)e*"-(™)V_a+p,„(t) 
-i3_(m)e-'^^-(™)>„(0a_a+}. 



(48) 



2 0.5 1 1.5 2 

FIG. 7: (Color online) Real and imaginary part of the coher- 
ence of the central spin for ao/ojo = 0.01 and A'' = 10. Blue FIG. 9: (Color online) The same as Fig. |8] for a longer inter- 
dots: numerical simulation. Green line: TCL approximation action time, 
according to Eq. (|49)) . 



where 




FIG. 8: (Color online) Population of the central spin for 
ao/ujo — 0.01 and A'^ — 10. Blue dots: Numerical simula- 
tion. Green line: TCL approximation according to Eq. (jSlf) . 



In the special case of uniform coupUngs (a^. = const) we 
have Ai = A2. Equation then reduces to the master 
equation derived in Ref. 17. 



The master equation (|l8|) can again be solved analyt- 
ically. The procedure is similar to the one outlined in 
Sec. ini We find the coherences, 



2N 



(49) 



, B_{m) 



N - 1 



{Al-Al)e 



^iQ-f {m)t 



+it 

and the populations, 



B+{m) B_{m) 



(50) 



iV+ 1 



N +1 



where 



8.42(A^+ 1) 



(1 - cos[n+{m)t]) . 



(51) 
(52) 



Involving only a sum over the quantum number m, the 
expressions (|49)) and ((5T|) can be evaluated numerically in 
a very efficient way. The comparison with the numerical 
simulation of the von Neumann equation demonstrates 
that also the TCL approach with the modified interac- 
tion picture yields a good agreement. The decoherence 
(see Fig. [7]) as well as the oscillations and the decay of 
the populations (see Fig. [8|) are very well reproduced by 
the simplified scheme. For longer interaction times the 
result ([51]) leads to revivals of the populations (see Fig. [9]) 
which arc due to the commensurability of the frequencies 
fi+(m). We stress that these revivals are neither present 
in the exact solution nor in the TCL approximation (|39p . 
Apart from these revivals the simplified approximation 
given by the master equation (j48p thus provides an accu- 
rate description of the decoherence and of the oscillating 
decay of the populations. 
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Finally we investigate the limit of a large number N 
of bath spins. For small values of the quantity (3 = 
2-\fN and large N the result (|5T|) can be approxi- 
mated by 



P+(i)« 1-/32 



1 



cos Wot 



(53) 



To obtain this expression one first expands the expo- 
nential in Eq. (|5T|) for small Ap°p and carries out the 
summation over m. Equation (|53p provides a good ap- 
proximation for a fixed /3 ^ 1 even for moderate A^- 
values. For example, the case of = 10 bath spins 
with ao/ciJo = 0.01 investigated above corresponds to 
/3 = 0.03. For this value of /3, we find that Eq. Q 
yields a good agreement for all N larger than about 10. 



V. CONCLUSIONS 

The appearance of a retarded memory kernel in the 
equations of motion, e. g. in the Nakajima-Zwanzig equa- 
tion, is often regarded as the characteristic feature of non- 
Markovian quantum processes. However, applications 
of the time-convolutionlcss projection operator technique 
show that strong non-Markovian behavior of open quan- 
tum systems can often be described by time-local master 
equations with an explicitly time-dependent generator. 
Although there is no theory of non-Markovian dynamics 
which allows a general assessment of the NZ and the TCL 
scheme, experience shows that in many physical systems 
the degree of accuracy achieved by both methods is of the 
same order of magnitude. In such cases the TCL method 
is clearly to be preferred because it only requires solving a 
time-local master equation. This does not mean that the 
TCL technique is always better than the NZ technique. 
The performance of these perturbation schemes strongly 
depends on the details of the system under investigation 



and on the chosen projection superoperator. There are 
examples of physical systems for which either the NZ or 
the TCL approach yields the exact result already in low- 
est order of perturbation theoryiii^. 

In the present paper wc have demonstrated the feasibil- 
ity of the TCL approach in the case of a nontrivial model 
describing a central spin coupled to a spin bath through 
nonuniform Heisenberg interaction. We have shown that 
the method indeed works in the strong field limit and 
provides a good approximation of the short and the long- 
time behavior of the coherences and the populations of 
the central spin. In addition, we have developed a TCL 
master equation that is based on a modified interaction 
picture and leads to a compact analytical solution for the 
central spin's density matrix which allows an efficient nu- 
merical computation. 

A possible approach to moderate or strong couplings 
is the analysis of the ordered cumulants of higher orders 
in the TCL expansion. However, it seems that a more ef- 
ficient strategy consists in the construction of more suit- 
able correlated projection superoperators. An example 
of this strategy is discussed in Ref. [13, where a certain 
correlated projection has been constructed for which the 
second order Nakajima-Zwanzig master equation leads 
to the exact dynamics of the population for the Hamil- 
tonian ([T|) with uniform couplings. It is of great rele- 
vance to develop possible extensions of this strategy to 
the nonuniform spin bath model analyzed here. 
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